Vortex fluctuations in underdoped Bi 2 Sr 2 CaCu 2 08+5 crystals 
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Vortex thermal fluctuations in heavily underdoped Bi2Sr2CaCu20g+a (T c =69.4 K) are studied using 
Josephson plasma resonance ( JPR) . From the zero-field data, we obtain the c-axis penetration depth 
Al,c(0) = 230 ± 10 /im and the anisotropy ratio j(T). The low plasma frequency allows us to study 
phase correlations over the whole vortex solid state, and to extract a wandering length r w of vortex 
pancakes. The temperature dependence of r w as well as its increase with dc magnetic field is 
explained by the renormalization of the vortex line tension by the fluctuations, suggesting that this 
softening is responsible for the dissociation of the vortices at the first order transition. 

PACS numbers: 74.60.Ec,74.40. Jg,74.60.Ge 



Vortex thermal fluctuations are considered essential in 
determining the (H, T) phase diagram of layered high 
temperature superconductors (HTS), and notably the 
first order transition (FOT) in which the ordered vortex 
crystal transforms to a liquid state without long range 
phase coherence || . Many scenarios, varying from vor- 
tex lattice melting described by a Lindemann criterion 
H to layer decoupling j|, g, g], all considering various 
degrees of coupling between the superconducting layers, 
have successfully been used to describe the position of 
the FOT in the (H, T)-plane. However, such fits to the 
FOT line have not been able to convincingly discrimi- 
nate between the different models. Here, we aim to do 
just this, through a direct measurement of the ampli- 
tude, as well as the field and temperature dependence of 
vortex thermal excursions in the vortex solid phase (or 
"Bragg-glass" 0) that lead to the FOT. 

For this study, we use the layered Bi 2 Sr 2 CaCu208+5 
(BSCCO) compound, in which vortex excursions can con- 
veniently be measured by the Josephson Plasma Res- 
onance (JPR) technique §, |, Briefly, the inter- 
layer Josephson current Jm can be measured through 
the JPR frequency u p i ~ Jm^ 2 , at which the equal- 
ity of charging and kinetic energy leads to a collective 
excitation of Cooper pairs across the layers. In turn, 
w 2 ;(B,T) = 07^(0, T)(cos(0„ in+ i)) intimately depends 
on the gauge-invariant phase difference <p n , n -\-i between 
adjacent layers n and n + 1 [n]. Here, (...) stands for 
thermal and disorder averaging. Thus, JPR is a probe 
of the inter layer phase coherence. The fluctuations of 
vortex lines created by a dc magnetic field applied per- 
pendicularly to the layers modify the relative phase dif- 
ference between adjacent layers and thus depress u) p i. In 
Bi2Sr2CaCu208+a, the ensemble of vortex lines should be 
described as stacks of two-dimensional pancake vortices. 
Thermal fluctuations shift the individual vortex pancakes 
with respect to their nearest neighbors in the c direction, 
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FIG. 1: Meandering of vortex pancakes along the vortex line 
in layered superconductors. Because of thermal fluctuations 
and disorder, pancakes (full circles) are shifted away from 
their equilibrium position (open circles). 



by a distance r„.„ + i = u n+ i — u„. Here u„ is the ab- 
plane displacement of the pancake vortex in layer n with 
respect to the equilibrium position of the stack it belongs 
to (Fig. 1). The wandering length of vortex lines, which 

can be then 
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is related directly to the JPR frequency 

defined as r w = (r^n+i) 1 2 |l3|, Below, we shall only 
consider temperatures above T — 42 K, at which vortex 
pinning (quenched disorder) is unimportant ^ p^[ . 

Underdoped BSCCO (T c = 69.4 K) single crystals were 
grown by the traveling solvent floating zone method in 
25 mbar O2 partial pressure at the FOM-ALMOS center, 
the Netherlands [jl5| . The samples were post-annealed for 
one week at 700° C in flowing N2. The advantage of using 
heavily underdoped BSCCO is that w pi (0,0) w 61 GHz 
turns out to be very low, which allows us to measure the 
vortex meandering over the entire vortex phase diagram. 
Samples A and B (cut from the same crystal) have di- 
mensions 1.35 x 1 x 0.04 mm 3 and 0.7 x 0.47x0.04 mm 3 , 
respectively. Another sample from the same batch was 
used to determine the temperature of the FOT (Fig. ||b, 
inset), by measuring the paramagnetic peak at the FOT 
with a miniature Hall probe magnetometer |]l6| . 

The JPR measurements were carried out using the cav- 
ity perturbation technique in the Laboratoire des Solides 
Irradies (on sample A) and the bolometric method in 
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FIG. 2: Field sweep experiment on sample A at T — 66 K in 
the TM012 mode of the cavity (/ = 22.9 GHz). The power 
absorbed (•, left axis) in the sample shows a maximum at 
Bjpr= 5.3 G at which ui = uj p i (open arrow). At the same 
field, the resonance frequency of the cavity (o, right axis) dis- 
plays a double-peak structure, indicated by the closed arrows, 
and a jump (arrow between dashed lines). 



the Institute for Solid State Physics at the University 
of Tokyo (on samples A and B). For the cavity pertur- 
bation technique, the sample was glued in the center of 
the top cover of a cylindrical Cu cavity used in the dif- 
ferent TM ii (i = 0, . . . , 4) modes. These provide the 
correct configuration of the microwave field at the sam- 
ple location, in which E r f\\c-axis and H r f w jL7|. The 
unloaded quality factor Qq is measured as function of 
temperature and field to obtain the power absorbed by 
the sample (Fig. 2). The bolometric method jOI con- 
sists in measuring the heating of the sample induced by 
the absorption of the incident microwave power when the 
JPR is excited |)|, [ijj. Furthermore, reversible magneti- 
zation measurements were carried out on sample A using 
a commercial superconducting quantum interference de- 
vice magnetometer in order to extract the London pene- 
tration depth \L,ab(T) for currents in the a6-plane ]l5| . 

Figure 3 shows the JPR frequency fjpR = Lu p i/2ir in 
zero field obtained by the above-mentioned methods on 
samples A and B. is proportional to the maximum 
interlayer Josephson current along the c-axis [pTj , 



(1) 

where Jm{B,T) — Jm (0, T)(cos(0 ni „ + i)) is the max- 
imum Josephson current, s = 1.5 nm is the interlayer 
spacing, e r the high-frequency relative dielectric con- 
stant and $0 the flux quantum. Using w p ;(0,T) = 
c/ 1 ^L,c{T)\fe r and e r =11.5 Q, we obtain the London 
penetration depth for currents along the c-axis, Al jC (T). 
When divided by the Ai ja h(T)-data from reversible mag- 
netization, this yields, without any model assumptions, 
the anisotropy parameter j(T) = Al, c (T)/ \L,ab(T), 
shown in the Inset to Fig. 0. Typically, at T = 0.5T C , 
Xl, c ~ 240 //m and \L,ab ~ 400 nm, so that 7 w 600, 
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FIG. 3: The JPR frequency in zero magnetic field for samples 
A (data below 40 GHz) and B (data above 40 GHz) vs the re- 
duced temperature T/T c . The bolometric method (•) and the 
cavity perturbation technique (o) have been used. We use a 
spline fit (solid line) in the extraction of the wandering length 
(see text). Inset: experimental temperature dependence of 7, 
obtained by the division of the experimental Xl,c(T) data by 
the A_L, a i,(T)-data from reversible magnetization. 

consistent with other data for the same material |j) . Note 
that 7 decreases as function of temperature. 

To analyze our JPR data in non-zero magnetic fields, 
we should divide ui p i(B,T) by the zero-field result de- 
picted in Fig. ||. In the absence of a model that sat- 
isfactorily describes uj p i(0,T) over the whole tempera- 
ture range, we resort to a spline fit to the experimental 
data. Next, we extract the vortex wandering length r w 
as follows. In the single vortex regime, at very low fields 
B < B.j = $ /A 2 7 , B < B\ = $o/4^A| jab , Bulaevskii 
and Koshelev derived fl3], 14 
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where the Josephson length Aj = 7s. We stress that this 
relation is meaningful only for small excursions r w < 

0. 6Aj, i.e. for (cos(0„,„ +1 )) = u 2 pl {B ,T) / u> 2 pl (0,T) < 

1. More generally, one expects an increase of 1 — 
(cos(<^ n>n _|_i)) with r w up to a plateau for large r w , 
as was found in recent simulations of the evolution of 
1 — (cos(0„. n+ i)) versus («)/<Zo ~ r w /ao for a pancake 
gas (ao = ^Qq/B is the intervortex spacing) ^0|. The 
numerical data show that 1 — (cos(</> rl . n +i)) is almost 
quadratic in r w for < 1 — (cos(</)„ i „ + i)) < 0.7 — 0.8, 
in agreement with Eq. (|2|), if the weak logarithmic de- 
pendence on Xj/r w is disregarded. Thus, we use 
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to obtain an approximation of the wandering length. 
Since r w = <(u„ +1 - u n ) 2 )V 2 = [2(u 2 - (u„ • u n+1 ))] 1 /2 ; 
one has, in the case of completely uncorrelated layers 
(e.g. for a pancake gas), r w — (2U 2 ) 1 / 2 = \[2u. Disre- 
garding the "anticorrelated" situation with u„-u n+ i < 0, 
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FIG. 4: 1 — (cos(0„,„+i)) vs temperature for different mag 
netic fields. We extract r w from these data using Eq.(S). 



correlations between pancake positions in different layers 
yield r w < y2u, i.e., r w is a lower limit for the root mean 
squared (RMS) displacement u of the vortex line. 

Figure 4 shows 1 - u 2 pl {B ,T) /uj 2 pl {0,T) = 1 - 
(cos(</> njn -|_i)) as function of temperature in different dc 
fields. The temperature dependence of the wandering 
length r w , obtained by applying Eq. (^), is represented 
in Fig. 5. For every field, we observe an increase of r w 
with T. At constant temperature, r w increases with mag- 
netic field, indicating a uj p i (-B)-dependence that deviates 
from the expected linear behavior of Eq. (||). Another 
interesting feature of the r w (T) curve is the break in 
the slope which appears at a field-dependent tempera- 
ture close to the FOT and above which all the r w curves 
merge into one. Alternatively, one may plot the same 
values of r w vs T/Tpor, where Tfot is the FOT tem- 
perature (Fig. 5b). Here, two regimes appear clearly, 
showing the correlation of the results with the FOT. 
For T < QMTfot, r w (T/T FOT ) roughly overlaps for 
all fields, whereas for T > 0.967>ot the curves deviate 
from each other. 

We now discuss the temperature and field dependence 
of r w in the vortex solid. The RMS thermal vortex dis- 
placement u can be obtained by equipartition of the as- 
sociated elastic energy with the thermal energy, U e i = 
C440q(u 2 /s) = fcflT. The vortex lattice tilt modulus 
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FIG. 5: (a) : Experimental r w vs T/T c for different mag- 
netic fields in strongly underdoped BSCCO. For B =27.5, 
22.4, 15.3, 12.4 and 9.9 G, arrows show the tempera- 
ture of the FOT. The dotted line shows the evolution of 
0.9s(fc B T7 2 /£os) 1/2 . Solid lines are fits to Eq. (|) with the 
single parameter a = 0.95, omitting the term 4/7r (a; 2 + 
for B > 15 G. (b) : r m vs T/Tfot- Solid lines are guides to 
the eye. Inset : phase diagram of a sample cut from the same 
crystal. Full and open triangles stand for the FOT and the 
irreversibility fields, respectively. 

the second term, introduced by the cutoff k max = 7t/r w , 
which corresponds to the smallest meaningful deforma- 
tion 21 ^2|. To proceed, we evaluate U e i at the typical 



vortex line deformation wavevectors parallel and perpen- 
dicular to the layers 
2-k/s. Writing Kq — 
equipartition yields 
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47T /ao, r 2 ^ = au 2 and x — ao/r w , 
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calculated by Koshelev and Vinokur |2lj and Goldin and 
Horowitz [ p2j | . consists of three terms: the nonlocal col- 
lective (lattice) term, the vortex line tension term, de- 
termined by Josephson coupling between layers, and a 
third term due to the electromagnetic dipole interaction 
between pancakes. Of particular interest here is the log- 
arithmic correction to the temperature dependence of 
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All parameters in Eq. (||), and notably e (T)/^/ 2 (T) — 
$o/47r/ioA 2 c , are known from experiment, which allows 
a direct comparison to the r w (T)-data. Very good agree- 
ment is obtained for the magnitude, the temperature, as 
well as the field dependence of r w for the lowest three 
fields, using the single free parameter a = 0.95. For 
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higher fields, Eq. (||) gives the correct magnitude of r w , 
but too weak a temperature dependence. However, ex- 
cellent fits of both the temperature and field dependence 
can be obtained for all fields, with the same a — 0.95, by 
omitting the nonlocal collective term, [i.e. 4/ir (x 2 + \) 
in Eq. (0)], see Fig. 5a. While the main temperature 
dependence of r w comes from the prefactor ~f 2 T/eo in 
Eq. ([)]) (dotted line in Fig. the behavior of r w in 
the vortex solid can only be understood as the result of 
the logarithmic correction arising from the softening of 
the line tension term by thermal fluctuations pl| , p2| . 
The field dependence, originating from Q z , explicitly in- 
dicates that vortex lines are correlated (line-like) on dis- 
tances that well exceed the layer spacing s. 

The experimental data can also be used to compare 
the terms entering Eq. (^). Deep inside the vortex solid, 
the line tension always dominates over the magnetic cou- 
pling and the nonlocal collective contribution. At very 
low fields (B < 10 G), the line tension term is largest all 
the way to the FOT. Eq. (|) then reduces to Eq. (40) of 
Ref. ^ with Q z ss irj /2ao instead of 2ir/s. At higher 
fields, the nonlocal collective contribution is expected 
to increase, eventually exceeding the Josephson coupling 
(line tension) term close to the FOT (for B > 20 G). 
Nevertheless, the very good fits obtained when only the 
line tension term is taken into account in Eq. (||) suggest 
that the line tension term always dominates C44 near the 
FOT. Moreover, we find that the electromagnetic cou- 
pling as well as the shear contribution to U e i are, under 
all circumstances, negligible. This renders Lindcmann- 
like j|] or dislocation-mediated (Kosterlitz-Thouless like) 
melting, as well as the vortex-line evaporation |6| scenar- 
ios very unlikely. Rather, the large thermal excursions of 
pancake vortices bring about the softening of the line ten- 
sion contribution to C44 for the large-wavevector modes 
that lead to the FOT. This would comply with recent 
measurements showing that vortex lattice order is not 
a prerequisite for the FOT |23]. For deformations with 
smaller wavevectors, Josephson coupling still contributes 
to the line tension even in the vortex liquid, leading to, 
e.g., the anisotropic vortex response to columnar defects 
in heavy-ion irradiated samples. 

Summarizing, we have carried out JPR measurements 
on heavily underdoped P^S^CaC^Os crystals. These 
data yield the c-axis penetration depth, the anisotropy 
parameter "f{T), and the wandering length of vortex lines 
r w . The observed temperature and field dependences of 
r w suggest that thermal fluctuations soften the Joseph- 
son coupling contribution to the tilt modulus for short 
wavelengths [^2| , which leads us to believe that this soft- 
ening drives the FOT. 
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